Effects of Inelastic Neutrino-Nucleus Scattering on Supernova Dynamics and 

Radiated Neutrino Spectra 
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Based on the shell model for Gamow- Teller and the Random Phase Approximation for forbidden 
transitions, we have calculated reaction rates for inelastic neutrino-nucleus scattering (INNS) under 
supernova (SN) conditions, assuming a matter composition given by Nuclear Statistical Equilibrium. 
The rates have been incorporated into state-of-the-art stellar core-collapse simulations with detailed 
energy-dependent neutrino transport. While no significant effect on the SN dynamics is observed, 
INNS increases the neutrino opacities noticeably and strongly reduces the high-energy tail of the 
neutrino spectrum emitted in the neutrino burst at shock breakout. Relatedly the expected event 
rates for the observation of such neutrinos by earthbound detectors are reduced by up to about 60%. 

PACS numbers: 97.60.Bw, 26.50.+X, 25.30.Pt, 95.55.Vj 



The simulation of core-collapse supernovae (SNe) is 
one of the great astrophysical challenges requiring so- 
phisticated computational skills and detailed input from 
various fields of physics (e.g., 0, Q)- It is wei l known 
that processes mediated by the weak interaction play an 
essential role for the collapse dynamics and the explosion 
mechanism [H, 13, [Bj]. However, inelastic neutral-current 
neutrino scattering on nuclei (INNS) , A + v — > A* + v' 
where A, A* describe the nucleus with mass number A 
before and after the scattering process, has not yet been 
satisfactorily considered in SN simulations. In this pro- 
cess energy is exchanged between neutrinos and matter, 
and hence it can potentially contribute to SN physics by 
(i) speeding up the thermalization of neutrinos with mat- 
ter after neutrino trapping during the collapse phase, (ii) 
changing the neutrino opacity which will in turn mod- 
ify the spectra of neutrinos released in the SN and to be 
observed by earthbound detectors, (iii) preheating the 
matter outside the shock front after bounce before ar- 
rival of the shock (iv) reviving the stalled shock, and 
(v) spallating nucleons from nuclei and thus modifying 
explosive nucleosynthesis. In this Letter we will report 
about SN simulations which consider INNS reactions and 
allow us to explore topics (i)-(iii) in detail. 

The effect of INNS in SN simulations has been inves- 
tigated previously in an exploratory study 0, approxi- 
mating the matter composition by a representative nu- 
cleus, 56 Fe. The reaction cross sections were based on a 
nuclear model for temperature T = 0, combining a trun- 
cated shell model evaluation of the allowed Gamow- Teller 
(GT) response to the cross section with estimates of for- 
bidden components derived from the Goldhaber- Teller 
model. The study concluded that INNS rates can com- 



pete with those of neutrino-electron scattering at moder- 
ate and high neutrino energies (E > 25 MeV), while they 
are significantly smaller for low E. No significant effects 
of INNS on the stalled shock by preheating the accreted 
matter were found in Q. 

Approximating the composition of SN matter by the 
ground state of the even-even nucleus 56 Fe is too sim- 
ple an assumption for the calculation of INNS interac- 
tion rates. At low and moderate neutrino energies, this 
neutral-current process is dominated by GT transitions 
mediated by the operator GTo = (ttq, where t±,Tq are 
the components of the isospin operator in spherical co- 
ordinates. The spin operator er can change the angular 
momentum of the initial state Ji by one unit and hence 
connects the 56 Fe ground state with spin/parity Ji = + 
only to final states in the same nucleus with Jf = 1 + . 
As the lowest 1 + state in 56 Fe is at an excitation en- 
ergy of E x — 3 MeV, there exists a threshold for inelastic 
neutrino scattering on 56 Fe (similarly on other even-even 
nuclei) and the cross sections are rather small for low 
neutrino energies. 

Supernova matter consists of a mixture of many nuclei 
with even and odd proton and neutron numbers. Since 
odd- A and odd-odd (i.e., with odd proton and neutron 
numbers) nuclei miss the strong pairing gap that lowers 
the ground state in even-even nuclei relative to excited 
states, and have usually Ji ^ 0, GT transitions from 
the ground state to levels at rather low excitation en- 
ergies are possible, reducing the threshold for inelastic 
neutrino scattering on the ground state and generally in- 
creasing the cross sections at low E. More importantly, 
SN matter has a non-zero temperature of order 1 MeV or 
higher, requiring the description of nuclei as a thermal 
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ensemble. This completely removes the energy threshold 
for inelastic neutrino scattering, because nuclei are with 
non-vanishing probability in excited states that can be 
connected to the ground state or final states at smaller 
excitation energies. Such scattering events correspond to 
de-excitation of the nucleus with the consequence that 
the final neutrino energy E' is larger than the initial en- 
ergy E. In fact, it has been demonstrated in Q that the 
consideration of T ^ effects can drastically increase the 
INNS cross sections for low neutrino energies. 

Our strategy to calculate the INNS cross sections ap- 
propriate for SN conditions is based on the fact that the 
matter composition at the relatively high temperatures 
(T > 1 MeV) is given by Nuclear Statistical Equilibrium 
(NSE). The computation of cross sections for individ- 
ual nuclei was described in where the allowed GT 
transitions for inelastic neutrino scattering on the ground 
states were calculated within the shell model. Due to the 
absence of experimental data for INNS this approach was 
validated by detailed comparison to precision Ml data 
from inelastic electron scattering, which for spherical nu- 
clei is dominated by transitions mediated by an operator 
proportional to GTo [1Q| . The contributions from forbid- 
den transitions, which become increasingly important at 
moderate and high E, were evaluated within the Ran- 
dom Phase Approximation (RPA). To incorporate the 
T ^ effects, our treatment distinguishes between 'up- 
scattering' (i.e. the final neutrino energy is smaller than 
the initial one; E' < E) and 'down-scattering' (£" > E) 
contributions to the cross sections. Up-scattering was 
treated approximately assuming the Brink hypothesis, 
i.e. assuming that the GTo and forbidden distributions 
on the excited states are the same as calculated for the 
ground state. The down-scattering contribution is ob- 
tained from the up-scattering contribution using detailed 
balance. 

In a SN environment the relevant INNS cross sec- 
tion d<j(E,E')/dE' is obtained by folding the cross 
sections d(ii{E , E') / dE' for individual nuclei with the 
appropriate abundance distributions, da(E,E')/dE' = 
J2i Yidai(E, E')/dE' = (Y\ Yi)(da(E, E')/dE'}, where 
the sum runs over all nuclei present and Yi denotes the 
abundance rii /rib of a given species (n^ and rib being the 
number densities of nuclei and baryons, respectively). 
Since the temperature is sufficiently high (T > 1 MeV) 
once the inelastic process becomes relevant, the nuclear 
composition is well approximated by NSE and, like in Q , 
we have calculated the abundances Yi from a Saha-likc 



NSE distribution, including Coulomb corrections [lll.ll2l|. 

Reference [9| has presented detailed cross sections for 
inelasic neutrino scattering on about 50 nuclei of the 
Z = 25-28 isotope chains. Additionally, using the 
method described above we have also included cross sec- 
tions for the Cr (Z — 24) isotopic chain. Importantly, 
for T > 1 MeV, all these calculations show little varia- 
tions between individual cross sections. Therefore, we 



have assumed that the average cross section over the full 
composition can be approximated by (dcr(E, E')/dE'} « 
Ei Yidai(E, E')/dE'/ J2[ Yi, where the sum is restricted 
to the pool of nuclei for which individual cross sec- 
tions have been calculated in We have determined 
a cross section table for a large variety of temperatures 
(0.517 < T[MeV] < 3.447), densities (10 s < p[gcm" 3 ] < 
6.31 x 10 12 ), and electron-to-baryon ratios (Y e , chosen 
density-dependent between 0.23 and 0.55), for a mesh 
of initial and final neutrino energies between zero and 
100 MeV. 
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FIG. 1: Normalized final-energy neutrino spectra for initial 
neutrino energies of 10, 30, and 50 MeV and T = 1 MeV, 
p = 6 x 10 11 g cm -3 and Y e = 0.397. The total com- 
position averaged cross section, (a(E)), for these energies 
is 0.48 x 10~ 42 , 38 x 10~ 42 , 243 x 10~ 42 cm 2 , respectively. 
For conditions of T = 2 MeV, p = 2.5 x 10 12 g cm" 3 and 
Y e = 0.275, the total composition averaged cross section is 



gies, showing a slight dependence on conditions. 

Figure [1] shows the normalized final-energy neutrino 
spectra at T = 1 MeV, p = 6 x 10 11 g cm" 3 and Y e = 
0.397 for three different initial energies. For E < 10 MeV 
down-scattering contributes significantly to the cross sec- 
tion; i.e., the de-excitation of thermally populated nu- 
clear levels produces neutrinos with E' > 10 MeV. Down- 
scattering becomes essentially irrelevant at higher neu- 
trino energies. For E = 30 MeV the cross section is 
dominated by the excitation of the GT centroid, giv- 
ing rise to a cross section peak around E' — 20 MeV. 
For even higher neutrino energies forbidden transitions 
contribute noticeably to the cross section. The peaks for 
E = 50 MeV correspond to excitations of the centroids 
of the GTo and dipole transition strengths. 

The SN calculations presented in this work were 
performed in spherical symmetry with the neutrino- 
hydrodynamics code Vertex (for details, see 13, [lj| )• 



The code module that integrates the nonrclativistic hy- 
drodynamics equations is a conservative, Eulerian imple- 
mentation of a Godunov-type scheme with higher-order 
spatial and temporal accuracy. The self-gravity of the 
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stellar gas is treated with an approximation to general 
relativity as discussed i n llal a nd tested against fully rel- 
ativistic simulations in [la , I IL 61 ] . The time-implicit trans- 
port routine solves the moment equations for neutrino 
number, energy, and momentum. It employs a vari- 
able Eddington closure factor that is obtained from it- 
crating to convergence a simplified Boltzmann equation 
coupled to the set of its moment equations. A state-of- 
the-art description of the interactions of neutrinos (v) 
and antineutrinos (v) of all flavors is included according 
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to 

We compare here the effects of INNS by simulatin g th e 
collapse of a 15 M Q progenitor star (model sl5a28 of [191]) 
with three different nuclear equations of state (EoSs) 
from Lattimer & Swesty (LS) [20|, Shen et al. [21(, and 
Hillebrandt & Wolff [22J, which are applied above some 
minimum density (~10 8 gcm~ 3 before shock breakout 
and 10 11 gcm -3 afterwards). At lower densities a mix- 
ture of ideal gases of e ± , photons, nucleons, and nuclei is 
used, and changes of the chemical composition are deter- 
mined by nuclear burning or through a 17-species NSE 
network. The high-density EoSs yield different time- 
variable abundances of neutrons, protons, a-particles, 
and a representative heavy nucleus, whose charge and 
mass numbers (Z, A) also differ between the three con- 
sidered cases. As a consequence, the INNS rates for 
v and v of all flavors are computed as the product of 
the tabulated pool-averaged differential cross sections, 
(da(E,E')/dE'), and an EoS-dependent abundance Ya 
of the representative heavy nucleus. Tests showed that 
insignificant changes of the SN results for a given EoS 
occur when the factor Ya in the rate calculation is re- 
placed by the abundance sum, £^1*, of all species of 
the NSE distribution. For reducing the dimensionality 
of the table, information of the cross section variation 
with the scattering angle was not stored. In the hydro- 
dynamic simulations we therefore made the approxima- 
tion that neutrinos colliding inelastically with nuclei are 
redistributed isotropically. 

The effects of INNS on the SN evolution and the prop- 
erties of the emitted neutrinos turn out to be very similar 
for all three employed nuclear EoSs. Despite causing a 
higher neutrino opacity by additional neutrino-nuclei in- 
teractions, INNS leads to a slightly stronger deleptoniza- 
tion and to a very small increase of the entropy in the 
homologously collapsing inner core. These effects are 
caused by up-scattering reactions in which high-energy 
neutrinos from electron captures lose part of their energy, 
thus producing additional heating of the stellar matter 
and escaping faster from the stellar interior because neu- 
trinos with lower energies possess a much smaller total 
interaction probability. The additional deleptonization, 
however, is so tiny that it reduces the collapse time to 
bounce and the enclosed mass of the shock formation 
radius only on a miniscule level. In contrast, due to 
the additional opacity, the number of electron neutri- 



nos radiated in the luminous burst that is released when 
the shock breaks out from the neutrino-opaque to the 
neutrino-transparent regime, is slightly smaller (on the 
third digit) when INNS is included. Moreover, these re- 
actions increase the neutrino-matter coupling and thus 
the total energy transfer rate in non-conservative scatter- 
ings (mainly on electrons and nuclei) ahead of the shock 
by typically a factor of 2-3 during the first < 50 ms af- 
ter core bounce (before the shock reaches a radius where 
the preshock density drops below the minimum value of 
our cross section table so that INNS was not taken into 
account any longer). Nevertheless, the preshock heating 
rates are large for too short a time to lead to consequences 
for the shock propagation and SN dynamics. The differ- 
ences with and without INNS remain smaller than the 
numerical resolution limit during all the simulated post- 
bounce evolution. 
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FIG. 2: Normalized v e number spectra radiated during the 
shock breakout burst as seen by a distant observer at rest. 
Results are shown for simulations with the three different nu- 
clear EoSs employed in this work. For better comparison 
of the strongly time-dependent spectra during this evolution 
phase, integration in a window of 8 ms around the peak lumi- 
nosity was performed. Inelastic neutrino scattering off nuclei 
(dashed lines) leads mostly to energy losses of high-energy 
neutrinos and thus reduces the high-energy tails of the spec- 
tra. The vertical lines mark the mean spectral energies. 

The most significant impact of INNS on the radiated 
neutrino spectra occurs in a time interval of about 15 ms 
around the shock breakout v e burst, when the preshock 
matter, which is composed mainly of heavy nuclei in 
NSE, has still a high density and therefore a fairly large 
optical depth for the escaping neutrinos. For this reason, 
high-energy neutrinos are efficiently degraded in energy 
space by frequent inelastic collisions with nuclei. While 
the mean spectral energy is reduced only by a modest 
amount (0.4-0.5 MeV), the high-energy tail of the emit- 
ted v e burst spectrum is strongly suppressed. The nor- 
malized spectra of electron neutrinos during a time inter- 
val of 8 ms around the maximum luminosity of the burst 
for simulations with and without INNS are displayed in 
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TABLE I: Electron neutrino cross sections for scattering off 
electrons and for charged-current interactions with nuclei of 
different detector materials, averaged over the SN spectrum 
of Fig. using the results obtained with the LS-EoS. The 
different cross sections were taken from: e [23 |. d [HI, 12 C [25^ . 
le O [2fJ, 40 Ar [27|, 56 Fe and 208 Pb @. 

Material (a) (10~ 42 cm 2 ) Reduction 





With INNS 


Without INNS 




e 


0.106 


0.110 


3% 


d 


4.92 


5.36 


8% 


12 C 


0.050 


0.080 


37% 


12 C (N ga ) 


0.046 


0.071 


35% 


16Q 


0.0053 


0.0128 


58% 


40 Ar 


13.4 


15.1 


11% 


56 Fe 


6.2 


7.5 


17% 


20Sp b 


103.3 


124.5 


17% 



Fig. [21 A similar effect of INNS is observed for P e and 
heavy- lepton neutrinos v x (see [l|) when their luminosi- 
ties begin to rise at the time the v e emission comes down 
from its peak but still remains clearly dominant for some 
ten milliseconds. 

Considering that the neutrino absorption cross sections 
on nuclei typically increase steeply with E, the reduc- 
tion of the high-energy spectral tail has important con- 
sequences for the detectability of the v c burst from SNe 
by experiments [lil • This is demonstrated in Table [J 
where we compare the relevant v c detection cross sec- 
tions calculated for the SN spectrum with and without 
INNS. In both cases the LS EoS was used; the other 
EoSs yield very similar results. The reduction depends, 
of course, strongly on the nuclear threshold energy. For 
12 C and 16 0, which are detector material in Borexino, 
MiniBooNe, SNO, and Super-Kamiokande, only neutri- 
nos with relatively high energies (E > 17MeV for 12 C 
and 15MeV for 16 0) can trigger charged-current reac- 
tions; consequently, the change in the SN spectrum re- 
duces the detection cross section by roughly 35% and 
nearly 60%, respectively. The data for 12 C(N gs ) include 
only the transition to the 12 N ground state. This is the 
only bound state in 12 N and the easiest transition to 
detect, because 12 N decays by e + emission after 11ms. 
The signal will be a e + e _ annihilation. We find a re- 
duction of 8% for the (v e ,e~) cross section on deuterons 
(the main SN detector material in SNO), 11% for 40 Ar 
(ICARUS), 17% for 56 Fe (Minos) and 208 Pb (OMNIS). 
In contrast, since the cross section for scattering off elec- 
trons increases linearly with E, its reduction is relatively 
small (only 3%), which is relevant for SNO and Super- 
Kamiokande. 

In summary, we have reported about the first SN simu- 
lations that include INNS, which allows for an additional 
mode of energy exchange between neutrinos and mat- 



ter. We found that this mode has little effect on the 
collapse dynamics and the shock propagation. However, 
INNS modifies the radiated neutrino spectra. In partic- 
ular, INNS strongly reduces the high-energy spectral tail 
of the v c burst at shock breakout. In turn, this notice- 
ably decreases the cross section for the observation of the 
burst neutrinos from future SNe by neutrino detectors. 

In Garching, this work was supported by DFG grants 
SFB/TR 27 and SFB 375. JMS acknowledges a grant 
from Fundagao para a Ciencia e Tecnologia. 
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